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r^ Abstract 

[~^ We consider the Curie- Weiss model at initial temperature < 13~^ < od in van- 

T-H ishing external field evolving under a Glauber spin-flip dynamics with temperature 

^^ < f3'~ < oo. We study the limiting conditional probabilities and their continuity 

r^ properties and discuss their set of points of discontinuity (bad points). We provide 

f*) I a complete analysis of the transition between Gibbsian and non-Gibbsian behavior 

as a function of time, extending earlier work for the case of independent spin-flip 

dynamics. 

For initial temperature /3~^ > 1 we prove that the time-evolved measure stays 
S Gibbs forever, for any (possibly low) temperature of the dynamics. 

In the regime of heating to low-temperatures from even lower temperatures, 

^vq < /3~^ < min{/3' , 1} we prove that the time-evolved measure is Gibbs initially 

J> and becomes non-Gibbs after a sharp transition time. We find this regime is further 

\l divided into a region where only symmetric bad configurations exist, and a region 

where this symmetry is broken. 

In the regime of further cooling from low-temperatures, /3'~ < /3~^ < 1 there 
is always symmetry-breaking in the set of bad configurations. These bad configura- 
tions are created by a new mechanism which is related to the occurrence of periodic 
^^ orbits for the vector field which describes the dynamics of Euler-Lagrange equations 

for the path large deviation functional for the order parameter. 

To our knowledge this is the first example of the rigorous study of non-Gibbsian 
phenomena related to cooling, albeit in a mean-field setup. 
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1 Introduction 

Non-Gibbsian measures are known to appear in many circumstances. Historically they 
were observed first in the context of position-space renormalization group transforma- 
tion and termed as so-called RG pathologies fllj. Later more and more examples were 
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discovered |2y4|8 , 19 , 24| which showed that the apphcation of many maps appHed to an 
infinite- volume Gibbs measure may result in similar "pathologies", meaning that the 
image measure is not a Gibbs measure anymore. When such a phenomenon appears 
it means that conditional probabilities of the image system will acquire long-range 
dependencies, at least for some non- removable configurations. Particularly interest- 
ing examples of infinite- volume transformations are coming from the study of dynam- 
ics [2| |14[jl6||23] . The first prototypical result in that direction is due to van Enter, 
Fernandez, den Hollander, Redig who considered an infinite-temperature (or high tem- 
perature) Glauber dynamics starting from an initial low-temperature Ising model on 
the two- or more-dimensional integer lattice. In particular they proved that a low- 
temperature initial measure in vanishing external magnetic field becomes non-Gibbs 
after sufficiently large times and stays non-Gibbs forever. This has to be contrasted 
with the simple fact that for independent dynamics, viewed on local observables, the 
time-evolved measure converges exponentially fast in time to the symmetric product 
measure. In fact such a phenomenon is possible since the Gibbs property (continuity 
property of conditional probabilities of the system) is to be tested in arbitrarily large 
volumes. Later more investigations for time-evolutions were performed. The general 
picture is that for very general dynamics and very general initial measures the time- 
evolved measures are again Gibbsian, for a sufficiently small time-interval |2|7| |17p9||20| . 
Long times however, even for simple dynamics offer the possibility for the emergence of 
non-Gibbsian measures. The discontinuities in the conditional probabilities which are 
responsible are produced by hidden phase transitions which pop up as a result of the 
conditioning procedure. Depending on the specific nature of the system there may be 
many mechanisms of such singularities [2, 14 . In this context continuous spin models 
are particularly interesting |6ll7,18 



While it is surprising that even the physically simple transformation of heating 
produces non-Gibbsian behavior it would even be more interesting to say something 
about cooling dynamics. More generally one would like to study a Gibbs measure fiQ 
for an initial Hamiltonian H which is subjected to a Glauber dynamics for another 
Hamiltonian H, which gives rise to a trajectory /ij where t denotes time. Glauber 
dynamics at low temeratures describes fast cooling or "quenching" . The question is to 
understand the behavior of fit, and in particular for which times it will be Gibbs. Since 
this is as yet too difficult on the lattice, we develop our results in mean-field. A mean- 
field system of Ising-type is called non-Gibbs if the single-site conditional probabilities 
depend in a discontinuous way on the magnetization of the conditioning spins [5 , 12 , 13 



17 . Investigations for mean- field models tend to reproduce the lattice results in many 
situations [17[[22] but often lead to an explicit knowledge of the parameter regions 
where Gibbsianness and non-Gibbsianness occur. Such an analysis has been performed 



for the Curie- Weiss model subjected to an independent spin-fiip dynamics in 14 . It 
was proved there that for initial high temperatures /3~^ > 1 the time-evolved measure 
is Gibbs forever, while for /3~^ < 1 there exists a sharp transition-time separating a 
Gibbsian from a non-Gibbsian regime. In the course of the analysis of that paper, also 
the phenomenon of symmetry-breaking in the set of bad configurations was observed 
which happens for the smaller range of initial temperatures below ,. In the present 
paper we build on that analysis but are able to extend the results to dependent spin- 
fiips according to a Glauber dynamics with an arbitrary other temperature P'~ . 

To understand discontinuous behavior of conditional probabilities for the time- 
evolved model at fixed time t one needs to look at the model resulting from the initial 



measure at time s = under application of the dynamics in the space-time region for 
times s between and t. The hidden phase transitions responsible for the non-Gibbsian 
behavior occur if there is a sensitive dependence of the model at time s = obtained 
from constraining the space-time measure to certain configurations at time s = t. If a 
small variation of such a constraining configuration leads to a jump in the constrained 
initial measure it will (generically) be a bad configuration for the conditional probabil- 
ities of the system at time t. Small variation means in the lattice case a perturbation 
in an annulus far away from the origin. Small variation means in the mean-field case 
a small variation of the magnetization as a real number. In the independent spin-flip 
lattice example of [4J the chessboard configuration was a bad one, correspondingly in 
the independent spin-fiip mean-field case of [14] the configurations with neutral magne- 
tization equal to zero were bad ones for large enough times. Moreover, configurations 
with non-zero magnetization also appeared as points of discontinuity for the limiting 
conditional probabilities, in a particular bounded region of the parameter space of initial 
temperature and time. This phenomenon was called biased non-Gibbsianness in p4j . 
The complete analysis for the mean-field independent spin-flip situation was possible 
since the constrained system on the first layer could be understood on the level of the 
magnetization. The relevant quantities could be computed in terms of the rate-function 
for a standard quenched disordered model, namely the Gurie- Weiss random-field Ising 
model with possibly non-symmetric random-field distribution of the quenched disorder. 

To deal with dependent-dynamics case a different route has to be taken since the 
dependence of the initial system on the conditioning is more intricate. As we will 
see, we will need to invoke the path large deviation principle for the dynamics with 
temperature f3'~ on the level of magnetizations. We will then have to minimize a cost 
functional of paths of magnetizations which is composed of the rate function along the 
path and an initial "punishment" term, which depends both on the initial Hamiltonian 
H and the dynamical Hamiltonian H, evaluated at the unknown initial point of the 
trajectory. The solution of the problem gives a surprising connection between path 
properties of the corresponding (integrable) dynamical system and Gibbs properties 
of a model of statistical mechanics. As a result we are providing a full description 
of the regions of Gibbsian and non-Gibbsian behavior as a function of time, initial 
temperature, and dynamical temperature. As a special case the previous results for 
infinite-temperature dynamics are reproduced (adding some geometrical insight about 
the behavior of typical paths) . Furthermore the solution reveals a new mechanism for 
the appearance of bad configurations in the region of cooling from low temperatures 
with even lower temperatures. These are related to periodic motion in the dynamical 
system. 

The present paper is to our knowledge the first one where Gibbs properties of a model 
subjected to a low-temperature dynamics are investigated and it will be challenging to 
see which parts of the behavior are occurring on the lattice. After the completion of our 
work we learnt about the preprint [3j where a large-deviation approach was proposed 
to understand dynamical transitions in the Gibbs properties for lattice systems, too. 
While there is a beautiful formalism available for path large deviations of empirical 
measures of lattice systems on an abstract level, explicit results are very hard, which 
underlines also the use of our present paper, and the necessity of future research. 

Moreover the questions and methods used should have interest also in models of 
population dynamics. In such models a population of N individuals, each individual 
carrying genes from a finite alphabet of possible types, performs a stochastic dynamics 



which can be described on the level of empirical distributions. Starting the dynamics 
from a known initial measure corresponds to an a-priori belief (prior distribution) over 
the distribution of types. Conditioning to a final configuration m' at time t corresponds 
to measuring the distribution of types. The occurrence of multiple histories leading 
to the same m' (which is responsible for non-Gibbsianness in the spin-model) has the 
interesting interpretation of a non-unique best estimator for the path explaining the 
present mix of genes. 



1.1 The model at time t=0 

We start at time s = with the Curie- Weiss Ising model in zero magnetic field 
at inverse temperature /3 whose finite-volume Gibbs measures on spin-configurations 



c^fi.Afl = i(^i)i=i,...,N G {-1, 1}^ are given by 



exp(^(^,^i(Ji)^^ 

^/3,Af(o-[l,Af]) = ^ (1) 

where the normalization factor Zi^^^- is the standard partition function. This model 
shows a phase transition at the critical temperature /3~^ = 1 in the limit where N — t- oo. 

1.2 The dynamics 

Given a configuration o-[i^n] = {(^i)i=i,...,N G {— 1) 1}^ of spins, we set Glauber dynamics 
on the level of the spins in such a way that it has the Curie- Weiss distribution with 
a (possibly) different temperature /3'~ as a reversible measure. We will call /3'~ the 
dynamical temperature. The generator of the system with A^ spins is given by 



N 
i=l j:jy^i 

where ah ^i denotes the configuration that is flipped at the site i 



^'l,N] 



''^■''^>r\ kS..^> '^' 



where we choose the rates to be 



J 



g±/3'm 



c(=p,m) = — . ,,o, ^ (4) 

cosh(p'?7ij — ?TT-smh(p'mj 

For fixed finite A^, we denote the corresponding time-evolved measure on {—1, 1}^ 
at time t, started from the equilibrium measure Hjs^n, by the symbol Hi3^/3'^t;N- It is clear 
that, for fixed A^, the time-evolved measure fJ^/3^i3'^t;N tends to the invariant measure 
under the dynamics, when 1 1 oo. 



1.3 The notion of Gibbsianness for mean-field models 

For a single-site spin ui G { — 1,+1} and a magnetization value for a system of size 
A^ — 1, that is m G {—1, — I + tv^T' • • • > ^~]v^' 1} '^^ consider the single-site conditional 
probabilities of the time-evolved measure in the volume A^ given by 

lp,P',t,N{(yi\fh) := ^/3,/3',i,Af(o-l|o-[2,Ar]), (5) 

where o"[2_Ar] is any spin-configuration such that fh = ]y^X]i=2'^j- -^y permutation 
invariance the right-hand side of (p]) does not depend on the choice of cr\2m ■ 

Definition 1.1 Let /3,f3',t be given. A point fh G (—1,1) is said to be good for the 
time-evolved mean-field model if and only if 

1. There exists a neighborhood of fh such that, for all a in this neighborhood the 
following holds. For all sequences a^ G {—1, —1 + ]v~l ' • • • > 1 ~ ]v^' ^} ^^^^ ^^^ 
property lim7v-|-oo oln = ol the limit 

7/3,/3',t(o-i|a) = lip 7/3,/3',i,Af(^i|aAf) (6) 

exists and is independent of the choice of the sequence. 

2. The function a i— t- 7/3,/3',t(o'i|a) is continuous at a = fh. 

Definition 1.2 The time-evolved mean-field model with parameters /3,P',t is called 
Gibbs iff it has no bad points. 

This definition has extensions to arbitrary local state spaces beyond finite types (see 
12,17]) where empirical magnetizations have to be replaced by empirical distributions 



in the definition. 

1.4 Main Theorem 

We are now in the position to give our main result. 

Theorem 1.3 Consider the time-evolved Curie-Weiss model with initial and dynamical 
temperatures I3~^, f3'~ . 

Then the following holds. 

1. Initial high temperature, any temperature of the dynamics. 

If P~^ ^ 1 then the time-evolved model is Cibbs for all t>0. 

2. Heating from an initial low-temperature, with a either high-temperature or a low- 
temperature dynamics. 

For any (5' there exists a value l^g^if^') < f3'~ (which is explicitly computable, see 
below) such that the following is true. 

Assume that < /3~^ < mm{l3'^^, 1}. 
(a) If Pskif^') < (3~^ then 



In 



1-/3 



the time-evolved model is Gibbs. 



• for aUO<t< t„Gs{^, /?') :- 4(13^ 

• for all t > tnGsif^y P') ^^6 model is not Gibbs and the time-evolved con- 
ditional probabilities are discontinuous at fh = and continuous at any 
fh ^ 0. 

(b) IfO< (3-^ < I3g^i(3') there exist sharp values < to(/3,/3') < ti(/3,/3') < 00 
such that 

• for all < t < to{P,f3') the time-evolved model is Gibbs, 

• for all to(/3, f3') < t < ti (/3, /3') there exists fhc = mc{(3, /?'; t) £ (0, 1) such 
that the limiting conditional probabilities are discontinuous at the points 
±fnc, <ind continuous otherwise, 

• for all t > ti(/3,/3') the limiting conditional probabilities are discontinu- 
ous at fh = Q and continuous at any in ^ 0. 

3. Gooling from initial low temperature. For f3'^ < (3^^ < 1 there exists a time- 
threshold tper{l3, /3') such that, 

• for all < t < tperi/3, P') the time-evolved model is Gibbs. 

• for all t > tperiP, /?') the model is not Gibbs and the time-evolved conditional 
probabilities are discontinuous at non-zero configurations fhc (and continuous 
at fh = 0). 
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Figure 1: Division between Gibbs and non-Gibbs area for low-temperature dynamics, 
the thick curve is obtained by computation, the dots are given by numerics 



Note that for high-temperature dynamics f3'^ > 1 the region 3 of initial temper- 
atures in Figure [T] is empty. Part 2 of the theorem generalizes the structure which 
we already know from the independent spin- flip dynamics /3' = (see [14] ) which is 
contained as a special case. This means that a symmetric (w.r.t. starting measure) 
bad point tuq = will appear after a sharp transition time if the initial temperature is 



not too low (see Subregion 2a). For lower temperatures (in Subregion 2b) symmetry- 
breaking in the set of bad configurations for the time-evolved measure appears in an 
intermediate time-interval: At the beginning of this interval a symmetric pair of bad 
configuration appears which merges at the end of the time interval. 

It is remarkable that the picture we observe in Region 2 is similar to the independent 
spin-flip case. This is even true for low temperatures /3' < 1 of the dynamics. As we 
will see, we can moreover compute the symmetry-breaking inverse temperature Psb in 
terms of /3' as the largest solution of the following cubic equation 

4/3iB + 12/3SB/3' - 6/3iB(l + /3') - /3'(3 + 3/3' - /3'') = (7) 

In the independent spin- flip case /3' = we get exactly /3~^ = |, which was already 
found in [14] . We will also give an explicit expression of the critical time in region 2a, 
for all /3'. 

In region 3 of cooling from an already low initial temperature we observe an entirely 
new mechanism for the production of non-Gibbsian points. These are related to periodic 
orbits of the flow of the /3'-dependent vector field which is created by the Euler-Lagrange 
equations obtained from the path-large-deviation principle for the given dynamics. 

1.5 Strategy of proof and phase-space picture 

To derive an expression for the time-evolved kernel 7/3,^' ,t it turns out that we need to 
look at path large deviations of the evolving empirical magnetization, on a fixed time- 
interval [0, t] with A^ as a large parameter, conditioned to end in a fixed magnetization 
m' S (—1,1). The path large deviation functional consists of two parts and can be 
viewed as a Lagrangian on the space of paths of magnetization on [0, i]. The first 
one is an integral over the time interval of a Lagrangian density depending on /3' as 
a parameter, and also on the magnetization variable and its time-derivative. Since 
the dynamics is started from an initial measure, the rate-functional in the LDP will 
contain also a second /3-dependent term "punishing" the choice of the (unknown) initial- 
condition. The solution of the corresponding path minimization problem will therefore 
depend on a balance between both terms. Such solution (or solutions, in case of multiple 
minima) will correspond to a most probable history path(s). Non- uniqueness of the 
solution makes possible a jump of the most probable history curve which ends at a 
prescribed final condition m' when one varies around particular values of m' . These 
particular values will become discontinuity points of 7/3,/3',t("i')- The problem of finding 
the most probable conditioned history path carries over analytically to the study of 
the evolution of a curve describing the allowed initial conditions for the magnetization 
and its velocity (depending on /?, /3') under the flow of the Euler-Lagrange equations 
(depending on /3'). Multiple histories show in this framework as multiple projections 
of the time-evolved curve in phase-space to the m-axis, and this will allow us to derive 
geometric insight as well as analytical and numerical results. As a warning we point 
out that the notion of "Hamiltonian" will always refer to a spin-Hamiltonian, not the 
Legendre transform of the discussed Lagrangian. 

The outline of the paper is as follows. Section 2 will be devoted to the derivation 
of the path large deviation principle, as well as the constrained large deviation prin- 
ciple involving the initial Hamiltonian and its relation to the time-evolved conditional 
probabilities. In Section 3 we discuss the solution of the variational problem in terms 
of the Euler-Lagrange equations giving rise to a time-evolved curve of allowed initial 



configurations. Section 4 provides more visual intuition for the system's behavior based 
on numerics. 
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2 Path large deviation principle and limiting conditional 
probabilities 

Before we start discussing a number of large-deviantions results it is appropriate to 
rewrite the finite-volume Gibbs measure (jTl) on spin-configurations crri^vi ^^ follows 



exp(-iVi7(mAf)) 



IJ'l3,N[Crii^N]) = ^— ; , (,oJ 



where rriN : (7[i^^] i— )• -^ J2i=i ^i i^ the function which sends a spin configuration to its 

empirical mean and H{x) := —2~ i^ ^^^ spin-Hamiltonian of the system. 

In this section we will provide an expression for the limiting conditional probabilities. 
This involves the large- A^ asymptotics for the paths of the empirical magnetization. Note 
first that, by permutation invariance, the continuous time process that is induced on 
the empirical magnetization is again a Markov chain. 

Namely, suppose that F : { — 1, — 1 + ;^, . . . , 1 — ^, 1} — )• M is a function on the 
possible magnetization values at size N and tjin is an empirical mean, then we have 
Li3\n{F o rriN) = {Lpi^j^F) o ttin with 



• N N \ 1 / / 2 

Lp,^NF{m) = ( — + — mjc(+,m- -j;^)\Py^ - J^] - F[m 





(9) 



How do typical paths for the unconstrained dynamics look for large A? Evaluating 
^E""o(F(mi)) = E""o((L^/,ArF)(mt)) for the observable F{m) = m for the expected 
value of the process started at ttt-q we have the identity 






(1 - mt)c{-,mt + —)-(! + mt)c{+, mt - —] 



(10) 



Taking the limit A — )• oo the magnetization concentrates on a deterministic path t i— )• 
m(t) which solves the ODE rh = (1 — m)c{—, m) — (1 + ?n,)c(+, m) or 

sinh(/3'?Ti-) — mcosh(/3'm) 
cosh(/3'?Ti) — msinh(/3'm) 



which has the largest solution of the mean-field equation m = tanh(/3'rn.) as a stable 
solution. In the case /3' = the equation reduces to the linear equation m(t) = —2m{t) 
which describes the relaxation of the magnetization to zero under the unconstrained 
infinite-temperature dynamics. 

Next we need to discuss a number of large deviation results which are needed to 
compute the limiting conditional probabilities. We begin as the first ingredient with 
the statement of the path large deviation principle for the dynamics with inverse tem- 
perature /?'. 

Theorem 2.1 Denote by Pp'^N the law of the paths (-ZAf(s))sG[o,tl of the magnetization 
for the continuous-time Markov-chain with generator Lpi^i\[. 

Then the measures P/s'^n satisfy a large deviation principle with rate N and rate 
function given by the Lagrange functional 



ip^Jpi{ip)= I jpi{ip{s),^p{s))ds 



with Lagrange density jpi (m, v) given by 



1 /e4/3'm(_l+^)2^2 + n +^)2^2_2e2/3'm(_l+^2W8_^^2) 

J3'{m,v) = -< 2 — 4 / K 



+ flog 
+ vlog 



■g-2/3'm /_;l ^ e2/''™(-l + m)-m 
4(-l + m) 



v + 



/e4/3'm(-_l + ^)2^2 + (1 + m)2^2 _ 2e2/3'm (_1 + m"^) (8 + v'^) 



For the special important case of non-interacting dynamics /3' = we write 
Jim,v) := JQ{m,v) 



- 2 - \/4 - 4m2 + 1)2 + V lo^ 



+ V4 - im? + 



2 -2m 



(12) 



(13) 



The proof of Theorem 2.1 will be sketched in the Appendix. This large deviation 
principle allows us to compute the large deviation asymptotics of finding the path of 
the magnetization jump process at finite N close to a given path (p{t). 

It allows us to compute the large deviation asymptotics of the probability to go from 
an initial configuration m to a final condition m' in time t by computing the value of 
the rate function in the minimizing path from m to m' . The minimizing path is found 
by solving the Euler-Lagrage equations with initial condition m and final condition m' . 



The second and more elementary ingredient we need is the static large deviation 
principle for the magnetization in the initial measure, the Curie- Weiss measure with 
inverse temperature /3. It reads as follows. 



Proposition 2.2 The distribution of magnetization m = ]v Si=i '^« w.r.t. the Curie- 
Weiss measure at inverse temperature /3 obeys a large deviation principle with rate N 
and rate function given by H{m) + I{m) where 

I{m) = — - — log(l + m)H — log(l - m) (14) 

is the rate function for the symmetric Bernoulli distribution. 

The Proposition is well known in the theory of mean- field systems. It follows from 
Varadhan's Lemma which states the following: Suppose a probability distribution sat- 
isfies a LDP principle with a known rate function and rate A^ and suppose we consider 
the probability distribution with density Ce~^^^^> relative to the first density. Then 
this probability distribution will satisfy a LDP with the same rate N and rate function 
obtained by adding H{m) to the first rate function and subtracting a constant. 

Combining these two ingredients we obtain the third statement governing the large 
deviation properties of the non-equilibrium system started in the inverse temperature 
P and driven with inverse temperature /3'. 

Theorem 2.3 Denote by Pp' ^p^N the law of the paths {zN{s))s^[o,t\ of the magnetiza- 
tion for the Markov- chain with inverse temperature f3' with initial condition distributed 
according to the Curie- Weiss measure fip'^N- 

Then the measures Pjs^/s^n satisfy a large deviation principle with rate N and rate 
function given by the Lagrange functional 

ip^H{^{0))+I{ipm + Jp,{^) (15) 

The knowledge of this compound rate function allows us to compute the large A^ 
asymptotics of the probability to find the system in a final magnetization m' at time t by 
computing the value in the rate function in the minimizing path to m' . Note that this 
time the optimization is also over the initial point m. This minimizing path ip is found 
by solving the Euler-Lagrange equations with final condition m' and an initial condition 
which is determined by another equation at the left-end point, which relates <^(0) and 
(p{0) in an /3- and /3'-dependent way, as we will see. We also call the corresponding 
curve in ((/^(O), (/j(0)) the curve of allowed initial configurations. 

Corollary 2.4 The conditional distribution of the initial magnetization m taken ac- 
cording to the law of the paths Pji^p^Nj conditioned to end in the final condition m' at 
time t, satisfies a large deviation principle with rate N and rate function given by 

Em'im,/3,P') = H{m) + Iim)+ inf Js'if) - Const{m') (16) 

<p:ip(0)—m, 

We are now ready to give our formula for the limiting conditional distributions of 
our model started at f3 and evolved with /?'. 



Theorem 2.5 Fix (3, (3' ,t,m' . Suppose the constrained variational problem (15) for 
paths if taken over the paths with fixed right endpoint ip{t) = m' has a unique minimizing 
path s I—)- m*{s;m' ,t). 



10 



Then the limiting probability kernels of the time-evolved measure fJ-i3^i3'^t;N have a 



well-defined infinite-volume limit ^j3 ^131 ^ti'liT^') i'n, the sense of (1.1) of the following form 



, , r, E.,=±ie'''^"'*^'-^"''^%t{ai,Vi;m',t) 

Here ps{cri,r]i;m' ,t) is the probability to go from ui G { — 1,1} at time s = to 
rji E { — 1, 1} at time s < t according to the Markov jump process on { — 1,1} which is 
defined by the time- dependent generator 

L{s; m', t)f{ai) = c(cri, m(s; m', t)){f{-ai) - f{ai)) (18) 

with rates which are obtained by substitution of the optimal path for the constrained 
problem for the empirical magnetization into the single-site flip rates. 

Proof. 

Take a sequence otv G { — !> —1 + j^zi, • • • , 1 — jfzi^ 1} with the property limaN = ol 

as A^ t 00. We denote by miv-i(s) = ]v^ Si=2 '^('^) *^^ empirical magnetization of 
the spins of site 2 to N . To prove that the promised form for the limiting conditional 
probabilities is correct we must show that 

moo P/3,/3',Ar(cJi(t) = r/'Jm7v_i(t) = Oat) 7/3,/3',t(??il«) 



where the r.h.s. is given by (17). 

Let us abbreviate the whole path (r7T,7v-i(s))o<s<t by the symbol x. Then, at finite 
A^, a double conditioning gives us the identity of the form 



■ P,P' ,N 



{ai{t) = 77i|mj^_j(t) = un) 



'^,3.^',iv(^l(*) = Vll"iiV-l(0 = oat) 



5-1 =±1 



the same with r/^ replacing r/i 

(20) 

We note that under our assumption on the solution of the constrained path large devi- 
ation principle the distribution ¥'p^pi^N{dx\mN^i{t) = a^) concentrates exponentially 
fast on the trajectory x* : s >-^ m*{s;m' ,t) as A^ tends to infinity. This collapses the 
outer expected value and simplifies the formula a lot. Next we have that whenever 
xn — )• X* we get 

„<Ti/9m*(0;rra',t) 

lim P^,^/,7v(cri(0) = ai\xN) = , .^ ^,„ ^ttt 

Ntoo 2cosh(pm*(0;m , tj) 

Finally we have that the single-site Markov chain describing the time-evolution of the 
spin at site 1, conditional on the path of the empirical mean of the other A^ — 1 spins 
and its initial value at time 0, converges to the Markov chain with deterministic but 



time-dependent generator (18). The corresponding transition probabilities converge to 



the limiting expression from the theorem and we have 

lim P/3,/3',Ar(o-i(t) =r/i|(Ti(0) =ai,XN) = pt{ai,rii;m' ,t) (21) 



This finishes the proof of (19). D 



11 



3 Phase-space geometry and multiple histories 

3.1 Euler-Lagrange equations and curve of allowed initial configura- 
tions 



Fix /3, /3', t, m' . We look at the constrained variational problem ( |15| ) taken over the paths 
(f with ip(t) = ml with the aim to find (the) minimizing path(s) s i— )• m*{s;m',t). It is 
known in the calculus of variations [10] (ch. 3, sect. 14) that a necessary condition for 
an extremum is given by the corresponding Euler-Lagrange equation and an additional 
free left-end condition of the form 

fsJvi^is), ip{s)) - j^{ip{s),^{s)) = for all s e [0, t] 

j^{if{s),^{s))-H^{^{s))-I^{^{s))\s=o = (22) 

ip{t) = m' 

where H denotes the initial Hamiltonian. Here we have dropped the subscript /3' for 
the function j{ip{s),ip{s)) and written subscripts to denote partial derivatives. It is 
straightforward to derive this set of equations by linear perturbation around the pre- 
sumed minimizing function ^{s) using a partial integration in the s integral. For more 
details, see the Appendix. 

The first equation is a second order ODE which has two unknown parameters which 
have to be determined by the second and third equation. We call the curve described 
by the second equation which gives a condition between initial point and initial slope 
of the solution curve the curve of the "allowed" initial configurations (ACC). We note 
that it is independent from the final value m' . 

Substituting the form oi j{ip{s),ip{s)), we get after a small computation the equa- 
tions 

((l+m)+e2/3'"'(l-m))^ 

g{m) 



m 



m 



s=0 

m[t) = m' 



(23) 



s=0 



with the function 



aim) = 2e2/3'-ii±II^I^i^!!!!:f2(l^^ 
^^ ' (i + ^) + e2™/3'(l-m) 



(24) 



describing the curve of allowed initial configurations. 

Here we have written m instead of (/^(s), and the dot denotes time derivative w.r.t. 
s. 

3.2 Typical paths for independent time-evolution 

Let us start with a discussion of the independent time-evolution. 
(i) For /3' = 0, /3 = 0, the system becomes 

rh{s) = Am{s) 

m{s) = 2m(s) (25) 

s=0 s=0 

m{t) = m' 

and the solution becomes m{s) = m'e'^^^~^'. This describes how a curve which is con- 
ditioned to end in m' away from zero is built up from the initial condition m'e~'^* close 
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to zero. 

(ii) For independent dynamics /?' = and initial inverse temperature /3 7^ the simpli- 
fied system is 



mis) 



m{s) 



s=0 

m{t) 



Am{s) 

g-2/3m(s)(i + m{s)) - e2/3™W(l - m{s)) 



s=0 



(26) 



m 



In this case the general solution is a linear combination of the e . Looking at the 
right-end condition one gets 



m 



(^s) = {m'-C2e^')e^^'-'^ + C2e^\ 



where C2 is a constant and must be determined by the left-end condition. This can be 
done numerically. 

It is possible to match the current approach with the one of [14J by plugging the 
solution curves with an initial condition ?tt,(0) = m* which are given by 



m[s) 



^ Of I 

m e — m 



pZiL p ZiL 



-2s 



+ 



m — m e 



*„-2t 



fZI-ZtJ £3 ZTj 



e^^sG [0,t] 



(27) 



into the rate function and carrying out the time integral explicitly. This gives 
E^i {m, (5, 0) = H{m) + I{m) 



+ - 4t + In 



+ ln 




+ 2 f m' In 



R - Cie-^* + Cae 



o2t 



1 — m' 



mln 



R-Ci + C; 

1 — m 



1-R- 2Cim'e-^* 1 + R- 2Cim 



l + R- 2Cim'e-2* I- R- 2Cim 
where R = ^/l - AC1C2, Ci 



2-f I 

me — m 



pZT, p ZiL 



, C2 



I —2/ 

m — me 



pZiL p ZiX, 



(28) 




Figure 2: Cost functional E^/^m, (3,0) and known function ^^,t,m'("^) ^^^ P' — ^1 ^'^d 
/3-1 « 1.744, i« 0.251 
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In the approach of |14j a related function cahed ^i3^t,m' (?7i) was obtained by Hubbard- 
Stratonovitch transformation, whose minimizers with a given conditioning (t, m') cor- 
respond to the most probable initial conditions. This provides an opportunity to check 
if the results of the present analysis done via path large deviations coincide with the 
approach employing the function ^i3^t,m'ii^)- 

It is known that the functions ^/3,t,m'("^) (3-2) and Em/im, /3,0) have the same set 
of extrema (see [21] in a more general context). In Figure^ is the plot of these functions 
(after normalization to have zero as a minimum) for the same set of parameters (/3, m' , t) 
which shows that the minima appear in fact at the same value. 

(ill) Let us next turn to the case of interacting dynamics /3' 7^ 0. In this case trajectories 
can only be obtained numerically. Before we go on, let us discuss in more detail the 
geometrical properties of the vector field and the allowed-configurations curve. 



3.3 Geometric interpretation of Euler-Lagrange vector-field and curve 
of allowed initial configurations 




Figure 3: Phase portrait with level curves and ACC, /3' 



Since the Euler-Lagrange density j{ip{s),(p{s)) (12) does not contain an explicit 
dependence on the time s, the generalized energy given by the Legendre transfrom of 



(12) is the system's first integral of motion 

j{(p{s), (p{s)) - <f{s)j^{ip{s), (p{s)) = C 



(29) 



This can be rewritten as 



g4/3'm(i_^)2^2 



+ (1 + mffn^ + 262/^'™ (1 - m^) (8 + th^) 



(1 + e2/3''"(l - m) + m)^ 
and explicitly solved for the velocity 



C 



(30) 



m=±JC+ 



16e2/3'r 



(m^ 



(1 -e2/3''»(m- 1) + m) 



(31) 
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Looking at the integral curves in phase space we get some geometric intuition. Let us go 



back to the notion of the ACC (24) on which all possible "allowed" starting conditions 



lie. In figure [3] there are several ACCs drawn which correspond to different values of 
/3, but the same value of the dynamical inverse temperature /3' = |, which is reltively 
low. The production of discontinuities of the limiting conditional probabilities will be 
related to the time-evolution of the curve of allowed initial configurations under the 
Euler-Lagrange vector field, as we will describe now. 

Let us first give a definition of a bad quadruple of initial temperature, dynamical 
temperature, time, and final magnetization in terms of dynamical-systems quantities. 
We start by defining candidate quadruples making use of the Euler-Lagrange flow in 
the following way. 

Definition 3.1 The quadruple {(3, (3' ,t,mpii) is called pre-bad iff there exists a pair 
"^0,1 7^ "10,2 of initial magnetizations s.t. the solution of the initial value problem of 
the Euler-Lagrange equations started in the corresponding points {mQ^i,g{mQ^i)) and 
(™'0,25 5(?Tio,2)) on the allowed- configurations curve for P , 13' has the same magnetization 
value niph at time t, that is 

m(f;mo,i,5(mo,i)) = m{t;mo,2,g{mo^2)) = mpb 

While this first definition refers only to the existence of overhangs of the time-evolved 
allowed-configurations curve, the next definition involves also the value of the cost (the 
large deviation functional together with the punishment term), which makes it much 
more restrictive. 

Definition 3.2 The pre-bad quadruple {l3, P',t,mi,ad) is called had iff the two different 
paths started at the corresponding m,o,i ^ mo,2 are both minimizers for the cost, i.e. 

E„,U^o,i,l3,P') = Em,^,{mo,2,P,f3') = miEm,Jm,P,(3') 



m 



We will exploit both definitions both to gain geometric insight as well as numerical re- 
sults. The important connection to non-Gibbsian behavior of the time-evolved measure 
lies in the fact that rribad of a bad quadruple will (generically) be a bad configuration for 
7;9,/3',t('|"^)- Indeed, to see this, let us go back to the explicit expression of the limiting 
conditional probabilities, given by 

113,13', t{Vl\m ) = ^ a.8m'(0-m't) ( ~ TTT ^^^1 

Note that the function m*(0;m',t) is not well defined for m! = m-bad itself since 
at time t there are two minimizing paths available, one from mo,i to m-bad a-nd one 
from mo,2 to m-bad- Varying however around mbad the paths will become unique and 
we might select the minimizing paths (and hence their initial points) by approaching 
the bad configuration from the right or left, obtaining (say) linim'irrn,^i'rn*{0]^ ,t) = 
mo,i and linim'^nn,^^ itt-* {0; m' , t) = mo,2- Note that we also expect that (generically) 
limm'imy,^APt{cri,fji;m',t) / limm"[nn,^iPticri,'ni;m' ,t). This follows since the pt are 
probabilities for two different single-particle Markov chains, one depending on the path 
starting from (?Tio,i,g(m.o,i))) the other one on the path starting from (mo,2, 5(^,0,2))- 
We note that, knowing the paths entering the pt^s, an explicit formula for pt in terms 
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of time-integrals can be written, and so, given (numerical) knowledge of the mini- 
mizing path, the 7/3,^',* (j?!!™.') can be obtained by simple integrations. Unless these 
two discontinuities compensate each other (which is generically not happening and 
which can be quickly checked by numerics) we will have that limm/j^^j^^^ 7/3,/3',t(^i|"i') 7^ 
liuim'^rrn,^^ 7^,/3',t(^i|"T'')- Consequently the model will be non-Gibbs at the time t. 

Conversely, if (/?, /3', t, rripb) is not bad, then m' 1— )■ 7^,/3',i(?7i|m.') is a continuity point. 
This follows since in that case all ?TT,'-dependent terms in (32) deform in a continuous 



way. So the absence of bad points (and a fortiori the absence of pre-bad points) implies 
Gibbsianness at (/3,/3',t). 



3.4 Time-evolved allowed initial configurations 

We just saw that non- Gibbsianness is produced by multiple histories which means in 
other words the production of overhangs in the time-evolved curve of allowed initial 
configurations. To get an intuition for this let us discuss the regions 2) and 3) of 
the main Theorem in more detail. Let us begin with the phase-space picture for the 
non-interacting dynamics /3' = 0. We are starting with the region 2a) of non-symmetry- 
breaking non-Gibbsianness i.e. | = /3^b(/3' = 0) < /3~^ < min{/3'~ ,!} = !. 




Figure 4: Non-symmetry-breaking mechanism, /3' = 0, /? ^ = 0.8 

The time-evolved allowed-configurations curve for t = tnGs(/5,/3' = 0) is shown at 
the left plot of Figure [4] where it acquires a vertical slope at zero. The right plot shows 
the time-evolved allowed-configurations curve for t > inGS (/?)/?' = 0) where it has two 
symmetric overhangs. In particular (/3, /3' = 0, t, m' = 0) is pre-bad. It is also bad, since 
the preimages of the upper and lower time-evolved allowed-configurations curve which 
intersect the vertical axis have paths with the same cost, by the symmetry of the model. 
Note that (/3,/3' = 0,t,m') is pre-bad for a whole interval of values of m' , but (as the 
study of the cost shows and as it was proved in (14] ) there are no other bad points. 
We note that m' = is easily checked to be indeed a bad configuration (discontinuity 
point) of 7^^^'=o,t("|^^') since there are no cancellations of discontinuities in this case, as 
we will explain now: Indeed, pt{ai,fji;m' ,t) does not depend on the trajectory of the 
empirical magnetization and is given by the independent spin-flip at the site 1 between 
plus and minus with rate 1, 



7p,P'=o,tivi\'m' 



2-^G\- 



:±r 



3O-i/9m*(0;m',t) 



Pt(o-i,??i) 



where Pi (+,-!-) 



i(l + e-2*), andpi(+, +) = pi(-, -) = l-pt(+, -) 



(33) 
l-Pi(-,+)- So, 
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a discontinuity under variation of m! is entering the formula only through r?T,*(0; m! , t), 
and hence m' i— >■ 7^,/3'=o,t(^i|^) is discontinuous iff m i— >■ m*{0;m' ,t) is discontinuous. 
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Figure 5: Symmetry-breaking mechanism, P' = 0,13 ^ = 0.4 
Let us now look at region 2b) of symmetry-breaking non-Gibbsianness i.e. /3" 



< 



/3s"b(/5' = 0) 

The left plot of Figure [5] shows the time-evolved allowed-configurations curve at 
t = to{f3, /3' = 0) where it acquires a vertical slope away from zero. The right plot shows 
the time-evolved allowed-configurations curve for to(/5, /3') < t < ti(/3,/3') where it has 
two symmetric overhangs away from zero. This means that (/3, f3' = 0, t, m') is pre-bad 
for a whole range of values of final magnetizations m' . Due to the lack of symmetry it is 
not clear to identify in the picture which of the (/3, /3' = 0, t, m')'s will be bad. It turns 
out that it is precisely one such value (/3,/3' = 0,t,mc), and this can be found looking 
numerically at the cost. 

Perturbations of these pictures stay true for /3'~ > 1, where they describe the 
only mechanism of non-Gibbsianness. Perturbations of these pictures also stay true 
for /3'~ < 1, but then there is also the Region 3 of the main theorem which de- 
scribes the cooling from an initial low temperature. We choose n = /3'~ < /3~^ = 
0.85 < 1. Then the vector field has periodic orbits which are intersected by the allowed- 
configurations curve, and the time-evolution will create overhangs and smear out the 
allowed-configurations curve over time. 




Figure 6: Non-Gibbsianness by periodicity, /3'' 



2^/3-1 = 0.85 



The left plot of Figure [6] shows the time-evolved allowed-configurations curve at 
t = tper(/3,/3') where it acquires a vertical slope away from zero inside the area of 
periodic motion. 
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The right plot shows the time-evolved allowed-configurations curve for a time t > 
to(/3,/3') where it has overhangs. Again, from the interval of pre-bad points, the bad 
point has to be selected by looking at the cost. When time gets larger more overhangs 
are created and the trajectory is smeared out. The corresponding potential function 
m I— 7- Efn'{m, /3, (3') will acquire more and more local extrema as t increases. Then, by 
finetuning of the m' while keeping the /3,/3',t fixed, equality of the depths of the two 
lowest minima can be achieved. Since the number of available minima is increasing with 
t we conjecture that there will be also an increasing number of bad m's which becomes 
dense as t increases. To prove this conjecture however, more investigation is needed. 

3.5 Emergence of bad points as a function of time 

The notion of a bad point can be viewed from two different standpoints. A pre-bad 
point in the time-space diagram is a point where two (or more) histories collide. If the 
costs computed along these paths are equal, then a pre-bad point is a bad point. In the 
phase space this means that the phase flow transported two (or more) points originally 
lying on the curve of allowed initial configurations to the same space-position within 
equal time but with different speeds. Two (or more) points have the same space-position 
if their projections to the m-axis are equal, as seen in Figures (|4|, ([s]), and (|6|. How 
can we identify analytically the first time t where time-evolved initial points from the 
curve of allowed initial configurations will obtain the same projection to the m-axis? As 
intuition suggests one has to look when the transported curve of allowed configurations 
aquires a vertical slope for the first time. This discussion brings us to the following 
computation. 

Writing u = m for the velocity, let us consider the flow m(t;mo,fo), v{t;mo,vo) of 
our system under the Euler-Lagrange equations, 

V = Jp'{rn) 

We take the curve of allowed initial configurations to be transported by the flow vq = 
9/3,13' (iT^o) where we write in short / = /^' and g = gp^p'- We are then interested in the 
projections to the m-axis of the time-evolved curves in phase space, that is the curves 
mo I—)- m{t;mQ,g{mQ)), as they evolve with t. Restricted to suitable neighboorhoods 
this curve becomes a function, and we view it as a potential function with state variable 
TTT-o and parameter t (keeping also f3,f3' as fixed parameters.) 

Doing so we see that the derivatives of the flow w.r.t. the initial conditions obey at 
the threshold time t that 

dm{t;mo,g{mQ)) dm{t;mo,vo) dm{t;mQ,vo) , 
U = r(3',i3[t,mo) := = 1 g [mo) 

(fm{t;mQ,g{mo)) 
{dniQY 

The first equation means that in the (m, v) plane the time-evolved curve will obtain a 
vertical slope which is clear by the interpretation of the variable mo as a parametrization 
of the curve of allowed initial configurations. 

Moreover we have that the second derivative will also vanish, since a minimum and 
a maximum of mo i— )• m{s] mQ,g{mo)) collide for s ^t, in a fold bifurcation. 
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3.6 The threshold time for non-symmetry-breaking non-Gibbsianness 
for dependent dynamics 

We can use these equations to obtain quantitative information about the threshold time 
for non-symmetry-breaking non-Gibbsianness also for dependent dynamics. For this it 
suffices to look at the dynamics locally around the origin {m, m) = (0,0) in phase space 
which is a stationary point for the dynamics independently of /?'. 
Linearizing f^'we get 

'J^) = (4(1-/3')' o)('^) ^^^^ 

Corrections are only of third order. The eigenvalues of the matrix are A1.2 = ±2(1-/3'), 
these eigenvalues are real and have different signs, so {m,rh) = (0, /^'(O)) = (0,0) is a 
saddle point. This ensures that the nature of solutions close to (0; 0) stays the same 
whatever /3' is taken. 

Let us now discuss the phase flow around the origin (0,0). At this point non- 
Gibbsianness without symmetry-breaking occurs, by the following argument. Suppose 
a symmetric pair of initial conditions (mo, v{mo)) and {—rriQ, t>(— m-o)) = {—mo, —v{mo)) 
is given which has the same time-evolved magnetization at time t. This corresponds 
to the fact that the transported curve will have overhangs at the points {0,vi{m)) and 
(0, —vi{m)). If we look at the phase portraits of the dynamics as a function of time we 
see that for times larger than but very close to the first time where this occurres the 
speed vi{m) will be very close to 0. It converges to when t approaches the transition 
time for Gibbsianness. Indeed the whole path was evolving in an arbitrarily small 
neighborhood of the origin and hence it suffices to look at the linearized dynamics. We 
also note that there is no need to look at the cost functional in this case, due to the 
symmetry of the paths. As time becomes larger than the transition-time (as in the right 
picture of Figure 4) the intersection points of the time-evolved curve with the vertical 
axis will move away from zero and so it would not be sufficient to use the linearization 
of the dynamics to compute the relation between bad magnetization values and time. 

Clearly the general solution of the linearized system is 

x{s) = Cie-^^'-^'> + C2e^^'-^'> (37) 

Putting the initial condition to be (mo, vq) the phase flow becomes 



m(s-ma Va) - 2(l-/30^O-^O ,~2(l-^')s 2(1 - /3')mo + T^O 2(1-/3'). 

m[s,mo,vo)- 4(1 _ ^/) ^ + 4(1-/3') ^ 

v{s;mo,vo) = !^C)^^a_^>o,-2(i-/3')s + vp + 2{1- ^')mo ^,^,_,,^, 



(38) 



Computing the function Fiii^p{t,mQ) (35) for this phase flow and setting it to zero, we 
solve it w.r.t. time t and get 

1 , 5'("io) - 2(1 - /3') 

In ^-rr-^. ^7^ ^ (39) 



4(1 - P') g'imo) + 2(1 - /3') 
Putting m-o = we obtain from this for the transition time 



' = 4(r^'"T^ <^"' 
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Figure 7: The upper branch shows the syninietry-breaking inverse temperature /3sb as 
a function of /3' 



By setting (3' 

iin(i-r') 



for the independent evolution in the last expression, the result t = 
is reproduced. We note that the transition time given by 



given m 
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formula (40) is positive only in the case when (3 > 1. 



To identify for which temperature- values the phenomenon of non-Gibbsianness with- 



out symmetry-breaking ends, let us look when the function (40) starts having several 



minima. In order to do this we compute the second derivative of (40) and put it equal 
to zero. This results in the equation 



4/3^ + 12/3/3' - 6/3^(1 + /3') - /3'(3 + 3/3' - /3'^) = 



(41) 



In the independent- dynamics case /3' = we get exactly /3 = |, which was already found 



in the paper 14 . 



3.7 Cooling and non-Gibbsianness by periodic orbits 

Let us specialize to the case of a low-temperature dynamics /3' > 1. In that case the 
phase space decomposes into the areas of periodic and non-periodic dynamics. The 
separatrix is given by (31) with C = 4. 



f±{m) 



±2 



(l + m)-e2/3''^(l-m) 
(1 + m) + e'^P'"^{l -m) 



(42) 



Note that the curve /+ (m) coincides with the curve of the "allowed" configurations ( 24 ) 
when /3' = /3. This means that it will be stable under the phase flow in that case. In 
particular the time-evolved curve will not acquire overhangs which corresponds to the 
fact that the time-evolved measure will be invariant under the dynamics and the model 
Gibbs. 

Note also that the negative branch of the separatrix coincides with the right-hand 



side of the ODE describing the unconstrained typical evolution (11) and so the inter 



section point with the ?TT--axis is given given by the biggest solution of the ordinary 
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mean- field equation m = tanh(/3'm). Let us first concentrate of the existence of pre-bad 
points, tliat is different initial points of the allowed-configurations curve leading to the 
same projection to the m-axis after time t. 

Now multiple overhangs are created if the allowed curve of initial configurations 
intersects the periodic motion area, as seen in Figure 6. Indeed this part of the curve 
will perform periodic motion and while doing so it will acquire more and more overhangs, 
filling out the part of the periodic motion area which is bounded by its extremal value 
of the integral of motion over time. It is now interesting to note for which temperatures 
this phenomenon can happen and this is the content of the following theorem. 

Theorem 3.3 (Non-Gibbsianness by periodicity) Suppose /S' > 1 and let m'Km'^) 
be the biggest solution of the mean-field equation for l3'{/3). Then the following is true. 

1. if 1 < /3 < 13' (or equivalently < 771-2 < ''^1) holds then 

• The curve of allowed initial configurations for /3, /?' has non-zero intersection 
with the (open) periodic motion area in phase phase for [5' . 

• Consequently there exists a threshold time tper{(3, /?') such that for all t > 
iper(/3,/3') there exists pre-bad {f3, 13' ,t,m')s. 

2. if 1 < 13 < 13' fails, there is either no periodic motion areas, or the curve of 
allowed-configurations has no intersection with them. 

Proof. Denote f = f^ (here we take the branch which bounds the periodic motion 
area from above), and the curve of the "allowed" configurations by g{x) (here x is used 
instead of m) so that we have 



g{x) = 2e 



{l + x) + e^^^'{l-xy 
,f, Al + x)-e^<P-P'\l-x) 
(l + a;)+e2^/^'(l-x) 



Previously it was mentioned that periodic motion arises only in the case /?' > 1, and so 
we will consider this along the proof, also w.l.g. we say that x > 0. Let us show what 
the condition 1 < /3 < /?' means and its equivalence to < mj^ < m^ First, we put 
f{x) = to determine the right border of the periodic motion area, and we get that it's 
given by the equation 

{l+x)~e^^'''{l~x) =0, 

which is equivalent to the mean-field equation for /3'. Second, consider f{x) = g{x) to 
determine their intersection point. This is simply 

(l + x)-e2^^(l-x) = 0, 

which is again the same mean- field equation, but for /?, where 7^2 has the same meaning 
as before. 

The allowed-configurations curve comes into the region of periodic motion and stays 
there when the following condition is satisfied 



fix) < g'{x) < fix) 



x=0 



x=Q 



x=0 
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Figure 8: Allowed-configurations curve for different /3 keeping /3' constant 

which turns out to be just equivalent to 

-(2/3' - 2) <2 - 4^3 + 2/3' < 2/3' - 2 



(44) 



or 1 < f3 < f3' . One can get an intuitive understanding of this mechanism from Figure 



4 Numerical results: Typical paths, bad configurations, 
multiple histories, forbidden regions 



Since the variational problem with fixed endpoint (23) can not be solved in closed form 



unless the dynamics is independent, let us now describe some of the key features which 
are seen in numerical study. 



For given conditioning (/3', /3, t,m') a solution of (23) with this set of parameters is 
called a history curve. Let us first discuss such curves for the example of the independent 
dynamics. Figure [9] shows on the right such history curves conditioned to end at time 
t at m' , for different values of m' . There is a jump in the optimal trajectory when we 
change m' = 0+ to m' = 0—. The associated cost functional at m' = 0, depicted on 



=0.8, t»l. 95608, m'=0 

o.o: 
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Figure 9: Symmetric forbidden regions 
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Figure 10: Non-symmetric forbidden region 
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Figure 11: Forbidden region for /3' = | 



the left, has two symmetric minima, and their minimizers are the two possible initial 
magnetization values. This is an example of a multiple history scenario. We call the 
regions showing on the right plot which cannot be visited by any integral curve forbidden 
regions. 



Figure 10 shows on the right history curves for the independent dynamics with 



a low initial temperature smaller than , where symmetry-breaking in the set of bad 
configurations takes place. We see on the right two discontinuity points m' and corre- 
spondingly two components of forbidden regions for the trajectory. The cost functional 
corresponding 

to the positive one of them is depicted on the right. Deformations of these pictures 
describe the phenomena for all temperatures of the dynamics, as long as the initial 
temperature is lower. 



Finally, Figure 11 displays history curves and cost functional at the critical condi- 
tioning for an example of cooling dynamics. 

Next, let us fix /3, /3' and describe the possible change of the set of bad configurations 
as a function of the time. Again we look at the independent dynamics first. 



The top line of Figure 12 has an initial temperature in which non-Gibbsian behavior 



without symmetry-breaking takes place. In the picture [l2b| we see the bad configurations 
m' as a function of the time s which were found numerically depicted by dots. Since 
m' = appears at a threshold time and stays to be the only bad configuration from 
that on, the graph of bad configurations is just a straight line starting at the threshold 
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Figure 12: Bad configurations as function of time (right) and initial points of trajectories 
(left) /?' = 



time. In the picture 12a we see the corresponding initial points of the history curves 
which are conditioned to end at m' . 



The lower line of 12 has an initial temperature for which non-Gibbsian behavior with 
symmetry-breaking takes place, in an intermediate time-interval. The right plot shows 
the corresponding non- negative branch of bad configurations m' . (By the symmetry of 
the model, taking the negative of these one obtains the full set of bad configurations.) 
The left plot shows the corresponding initial points of the history curves which are 
conditioned to end at the non-negative bad configurations m' on the right. 



Finally, figure 13 displays the time-evolution of bad configurations and their initial 
points for a low-temperature dynamics. The lowest line corresponds to heating from 
very low initial temperature and shows non-Gibbsianness with symmetry-breaking at an 
intermediate time-interval. The middle line corresponds to heating from an intermediate 
lower temperature and shows non-Gibbsianness without symmetry-breaking. These 
two mechanisms are known from high-temperature dynamics. Figures 13a and 13b 
correspond to cooling and shows data from the region of periodic orbits. 

Applying numerical integration of the Euler-Lagrange equations from initial condi- 
tions chosen on the allowed-configurations curve, check for intersecting trajectories and 
numerical computation of the cost function we can get (numerical approximations to) the 
array of bad quadruples (/3, /3', t,mpb), augmented by the possible initial points. With 
this procedure we rederived the Gibbs-non-Gibbs phase diagram for /3' = (which was 
obtained earlier in [l4|). Based on it we can draw the Gibbs-non-Gibbs phase-diagram 
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Figure 13: Bad configurations as function of time (right) and initial points of trajectories 
(left) - low-temperature dynamics /?' = § 

at any dynamical temperature f3' . An example for this was presented in the Introduction 
of the present paper in figure [T] for a fixed low dynamical temperature. 
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5 Appendix 

5.1 Sketch of proof of unconstrained path large deviation principle 

Let us first consider first a simpler Markov jump process for tlie magnetization witli 
transition rates wliich to not depend on the state m of the process. The corresponding 
path large deviation principle can be built up as follows. The first ingredient is the large 
deviation principle for the Compound Poisson process. 

Proposition 5.1 Denote by Rt a Poisson process with rate 1. Denote by ^j, i = 1, 2, . . . 
a sequence of i.i.d. random variables with exponential moment generating function 
h{X) =Ee^^K 

Denote by Z{t) = X]j=i d ^he associated Compound Poisson process. Define the 
rescaled paths by Zj^lt) = j^Z{Nt). 

Define 

J{v) = sup(fA - h{\)) + 1 ^45-) 

Then, at fixed t, as N "[ 00, the distribution of the variable Zj^jit) satisfies a LDP in M 
with rate function tJ {j) and rate N . 

This is a known theorem ilj but it is instructive to see the proof for the sake of com- 
pleteness. 

Proof: Let us look at the logarithmic moment generating function Ajv(A), defined 
by 



^ ^ to ^ ^ ^' (46) 

= e(-^*+^*K^)) 

Recall the Gartner-Ellis theorem (Theorem 2.3.6. in [Ij, page 44) which states the 
following. Assume that limAr -^Ajv(A^A) =: A(A) exists. Then the distribution of the 
variable Zfyf{t) satisfies a LDP on M with rate function A*{v), which is the Fenchel- 
Legendre transform of A(A). In our case we have equality even at finite A^ of the form 
A(A) = -^A7v(A^A) = t{h{X) — 1) and the Legendre transform gives 

A*{v) = sup(«A - t{h{X) - 1)) = tJ(^) (47) 

A * 

So the theorem follows. D 

In the next step we go from one-dimensional large deviations to large deviations of 
finite-imensional marginals in the path space of the Compound Poisson process. This 
way of arguing corresponds to [I] Lemma 5.1.8, page 178, which is a step to prove 
Mogulskii's theorem. Recall that Mogulskii's theorem states that the paths of empirical 
averages of the form -^ Yli=i ^i satisfy a LDP. 

Proposition 5.2 For any decomposition B = {0 < ti < t2 < ■ ■ ■ < t^ < t} of the 

time-interval [0, t] and path f : [0, t] — t- M denote by irsf the projection 

7^Bf={fih)J{t2),...,fitk)) (48) 
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Then the corresponding image measures ttb{Pn) = Pn ° t^b ^f ^^^ rescaled Compound 
Poisson process satisfy a large deviation principle in M with rate N and rate function 

Jb{z,, ...,z,) = ti j(^) + (t2 - h)j{'^) + • • • + (tfc - tfc-i)j( f ~f-M 

\ti/ \t2-ti/ ^tk - tk-i y 

(49) 



where J is defined in (45). 



The proof follows from putting together the result for the one-dimensional distributions 
as above. From here one gets the analogue of Mogulskii's theorem (see llj Theorem 
5.1.2, page 176). 

Theorem 5.3 Denote by P/v the law of the rescaled paths by Z]\j{s) = jjZ{Ns) of the 
Compound Poisson process as above, for < s < t. 

Then the measures Pn satisfy a large deviation principle with rate N and rate func- 
tion given by the Lagrange functional C{ip) = L J{ip{s))ds where J is defined in (45). 



We do not give a full proof here, but note that it follows by taking a supremum 
over the finite decompositions, and invoking the Dawson-Gartner theorem, as explained 
in [I]. 

Up to this moment we have only treated the Compound Poisson process with con- 
stant rates, which is not the case here, since we consider state-dependent spin-flip dy- 
namics, meaning state-dependent rates. We need to justify rigorously that we can 
replace the contribution to the integral over the Lagrangian density for the infinitesimal 
time-interval ds of the form j{(p{s))ds by a term jpi[ip{s), (p{s))ds if the transition rates 
of the Markov chain depend on the state '^{s). In order to do that a comparison result 
is needed which compares a Markov chain with constant rates with the original Markov 
chain with state-dependent (but bounded) rates on the level of the logarithm of the 
exponential moment generating function, on small tiine-intervals. We are grateful to 
Prank Redig for communicating the following result to us which is an essential ingredi- 
ent. Informally, the following lemma states that a jump process z^it) on the discrete 
space with constant jump-up, jump-down rates c± does not differ much from a process 
mN{t) with state-dependent rates c±{mN{t)) when the time interval [0, At], where both 
processes are considered, is sufficiently small and N is sufficiently large. This is due to 
the fact that the state of mN{At) cannot change a lot if At is small and N is large. At 
small times mi\i{t) can make not much more jumps than A^At jumps of a small height, 
therefore the state-dependent rates c±{mN{0)) and c±{mj\f{At)) will not vary much. 

Lemma 5.4 F. Redig's useful lemma. Denote by ZM{t) the Markov process on the 
discrete space { — 1,-1 -|- ;^,...,1 — ;^,1} started at mo with constant non-zero rates 
c±(77io) to go up of down by one step, and by mN{t) the true process started at the same 
point mo with state- dependent rates c±{m) given by ([9]). Then 

,. 1 ,. 1 , EexpJVAzjv(t) 

lim - sup lim sup — : log — — ^-^ = L) ( 5U ) 

m t m Ntoo N Emex.pNXmN[t) 

The full proof will appear elsewhere, along with generalizations to more general local 
state spaces. Employing the lemma and going through suprema over finite partitions 



again, the proof of Theorem 2.1 is obtained where we still need to identify the form 
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of the Lagrangian density. Let us start again with a Compound Poisson process with 
jumps of size 2 with the distribution P{^i = 2) = p, P{(,i = —2) = 1 — p where p is fixed 
in the beginning. If we denote again /i(A) = Ee''^^^ and J{v) = sup;^(vA — h{X)) + 1 we 
have that 

J» = ^ln(^^±^ ^^-^^^' + ^y ^Vl6p-16p^ + .^ + l (51) 



as a solution of a quadratic equation shows. Choosing the rates in the rescaled Com- 

W (or down by -^; 



pound Poisson process to go up by ^ (or down by — ^) to match the rates in the 



generator L^/^^ we are led to choose 

P^'i^) = 2B'mn N I n I A (^^'' 

and this explains the form of the Lagrangian density Jj3'{m,v) = Jp ,(rn){v), after a 
small computation. This concludes our treatment of the proof. D 

5.2 Free end-condition 

To obtain the necessary condition ( |22| ) for an extremum of the variational problem 

ip^H{^{0))+I{ip{0)) + Jp,{^) (53) 

with ip{t) = m' , use the standard procedure in calculus of variations adapted to the 
problem with a free left end: Consider a perturbation ip{s) + eAip{s) around the ex- 
tremum ip{s), with a function Ac/?(s) obeying the constaint Aip{t) = at the end-point 



but no constraint on Ac^(O) at the initial point. Plug ^{s) + eAip{s) into (53), expand 
to linear order in e, and demand that the terms proportional to e vanish. Using partial 
integration under the s-integral one arrives at the Euler-Lagrange equation for s in the 
interval between and t, and the additional free end condition at the initial point, the 
latter one following by demanding that the terms proportional to Aip{0) have to vanish. 
Alternatively, the problem can be reformulated in terms of a problem with different 
Lagrange density but without initial punishment term, via incorporation of the ini- 



tial term into the integrand. From here, we refer to [10 where free-end problems are 
discussed. 

5.3 Hamiltonian and Lagrangian formalism 

Following a suggestion of a referee, let us remark that, as an alternative to the derivation 
in terms of approximations by compound Poisson processes, the form of the Lagrangian 
can also be obtained going through the formalism of \9\, see Chapter 1.4. Let us briefly 
sketch this procedure for the convenience of the reader. The jump process for the 
magnetization has a generator 

/ 2 2 \ 

AN9{m) = N[pp'{m){g{m+ — ) - gim)) +(1 - p[si{m)){g{m - — ) - g{m))j 

ff z \ (^^) 

=:N j yg{m + ^) " 9{m))v{'m', dz) 
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with ppi{rn) given by (52). This generator is of the form treated in p\ with the obvious 



identification of rj{m,dz). From here one defines an operator 1-L by the corresponding 
action on functions / of the magnetization of the form 

{nf){m)= /(e^'(™)^-l)r?(m,<iz) 

J (55) 

= pp,{m){e^f'^'^^ - 1) + (1 -p;3,(m))(e-2/'("^) - 1) 

Following the formalism and replacing f'{m) by a momentum variable A one defines the 
corresponding governing Hamiltonian function l-L{m,X) (or generalized energy) of the 
dynamics (not to be confused with the spin-Hamiltonian) 

H{m,\) = pp>{m){e^^ - 1) + {I - pp,{m)){e-'^^ - 1) (56) 

Performing a Legendre transform we arrive at the Lagrangian density 

jfji{m,m) = sup(Arh — 'H{m, A)) 

A 

which governs the path large deviations (as stated in [9], page 12). This procedure 
was also employed for explicit computations in the infinite-temperature case in [3]. 
It is a computational exercise to verify that this approach reproduces the form of the 
Lagrangian density previously given. Note also that the integral of motion we introduced 



in (29) is identical to the generalized energy 7i{(f{s), (p{s)). 

Let us remark that the study of the corresponding Hamiltonian-Jacobi equations is 
equivalent to the study of the Euler-Lagrange equations. 
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